clc;
clear all;
close all;
% Define the differential equation
syms y(x);
eqn=input('enter the differential equation')
cond=input('Enter the Boundary Condition')
yg= simplify( dsolve(eqn))
yp = simplify(dsolve(eqn,cond))
ezplot(yp)